function M = printResults(results,reb,mpow,specif)

beta = {'beta:const', 'beta:terrain', 'beta:gdp','beta:dem', 'beta:pop'};
phi1 = {'phi1:US', 'phi1:UK', 'phi1:FR', 'phi1:RU', 'phi1:CH', ...
    'phi1:terrain', 'phi1:gdp','phi1:dem', 'phi1:pop'};
phi2 = {'phi2:US', 'phi2:UK', 'phi2:FR', 'phi2:RU', 'phi2:CH', ...
    'phi2:terrain', 'phi2:gdp','phi2:dem', 'phi2:pop'};
gamma1 = {'gamma1:US', 'gamma1:UK', 'gamma1:FR', 'gamma1:RU', 'gamma1:CH', 'gamma1:dist', 'gamma1:allies', ...
    'gamma1:colony', 'gamma1:war'};
gamma2 = {'gamma2:US', 'gamma2:UK', 'gamma2:FR', 'gamma2:RU', 'gamma2:CH', 'gamma2:dist', 'gamma2:allies', ...
    'gamma2:colony', 'gamma2:war'};
chi1 = {'chi1:dist', 'chi1:allies', 'chi1:colony', 'chi1:war'};
chi2 = {'chi2:dist', 'chi2:allies', 'chi2:colony', 'chi2:war'};
delta = {'dS:US-UK', 'dS:US-FR', 'dS:US-RU', 'dS:US-CH', 'dS:UK-FR', 'dS:UK-RU', 'dS:UK-CH', 'dS:FR-RU', 'dS:FR-CH', 'dS:RU-CH', ...
    'dO:US-UK', 'dO:US-FR', 'dO:US-RU', 'dO:US-CH', 'dO:UK-FR', 'dO:UK-RU', 'dO:UK-CH', 'dO:FR-RU', 'dO:FR-CH', 'dO:RU-CH'};
lambda = {'l:Reb_ord', 'l:US_ord','l:UK_ord', 'l:FR_ord', 'l:RU_ord', 'l:CH_ord', ...
    'l:Reb_binary', 'l:US_binary','l:UK_binary', 'l:FR_binary', 'l:RU_binary', 'l:CH_binary', ...
    'l:ExpInterv', 'l:ExpIntervSG','l:PureStrategy', 'l:War', 'l:ParetoDominated'};

rowNames = [beta(reb{1}),phi1([1:4,4+mpow{1}]),phi2([1:4,4+mpow{2}]),gamma1(1:5),gamma1(5+reb{2}),...
    gamma2(1:5),gamma2(5+reb{3}),chi1(mpow{3}),chi2(mpow{4}),delta,lambda(specif)];

varNames = {'Estimate', 'StdError', 'pValue'};

M = array2table(results,'RowNames',rowNames,'VariableNames',varNames(1:size(results,2)));
